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ABSTRACT 


Two different versions of the Green’s function for the 
scalar wave equation In weakly curved spacetime (one due to DeWltt 
and DeWltt, the other to Thome and Kov4cs) are compered and con- 
trasted; and their mathematical equivalence is demonstrated. Then 
the DeWitt-DeWitt Green’s function is used to construct several 
alternative versions of the Thome -Kovacs post-linear formalism for 
gravitational-wave, generation. Finally it is shown that, in calcu- 
lations of gravitational bremsst rah lung radiation, some of our ver- 
sions of the post-linear formalism allow one to treat the interact- 
ing bodies as point masses, while others do not. 



I. INTRODUCTION 


In the first paper in this series Thome and Kovacs (1975) [cited 
henceforth as TK] developed a "plug-in-and-grind*,* "post-linear” formalism 
for calculating the gravitational radiation from accelerated systems with 
weak internal gravitational fields but arbitrarily large internal veloci- 
ties. Their method is applicable to a class of problems which heretofore 
have not been amenable to analysis. The application to the gravitational 
bremsstrahlmg problem is of particular importance. Central to the devel- 
opment of their formalism is the determination of an appropriate approxi- 
mation to the well-known exact Green's function for the scalar wave equa- 
tion in curved spacetime developed by DeWitt and Brehme (1960) . 

The purpose of this note is threefoldt First , to show that an 
approximate form of the exact Green's function developed previously by 
DeWitt and DeWitt (1964) [cited henceforth as DD] has exactly the same 
mathematical content (within the constraints of the approximations ixsed) 
as that constructed by TK, although the two representations differ signif- 
icantly in mathematical form and physical interT>retation, and have slightly 
different realms of validity. Second , to use the Green's function of DD 
in the formulas of TK to establish several alternative representations of 
the TK post-linear gravitational-wave-generation formalism. Third , to 
determine, in the case of the gravitational bremsstrahlung problem, which 
representations of the formalism permit one to treat the interacting ob- 
jects as point masses rather than as extended bodies, and which do not. 

In §II wc review and contrast the DD method of approximation with 
that of TK. The equivalence of the two approximations is demonstrated in 
§III. In §IV the alternative representations of the "plug-ln-and-grlnd" 
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formalism are given. Section V discusses the polnt-*mass approximation to 
the bremsstrahlung problem. 

In the following we shall deal with three sets of 2-tensors, namely, 
those associated with the event pairs (x,x'), (x,x") and (x",x'). To 
avoid confusion we shall distinguish the various 2-tensors by bars and 
tildes as follows (cf. fig. 1): 

A - A(x,x'), A - A(x,x") and A - A(x",x’) . (1) 

Further, we shall let the Indices on 2-tensors designate the events to 

u u’ 

which they refer (e.g,: A = 9A/9x^ , A , = 9A/9x etc.). 

Finally, we shall use geometrlzed units (c - G •* 1) . 


II. WEAK-FIELD GREEN'S FUNCTIONS 
a) Foundations for the Analysis 

The defining equation for the Green's fxjnctlon, corresponding to the 
scalar wave equation in curved spacetime, is given by 

g'*'’ , (2) 

where g is the determinant of the metric tensor g^^ , semicolons indi- 

4 

cate covariant differentiation, and 5 (x,x') is the four-dimensional bl- 
denslty Dirac delta function 


6^ = 6^(x,x') 


6(x° - 


X ) 


6(x^- 




6(x^- 


X ) 


6(x^- 


X ) 


6^(x' ,x) . 
(3) 


In the absence of the crossing of light-cone geodesics, the exact solution of 
equation (2) (due to DeWitt and Brehme [I960]) consists of a "direct part" 
and a "tail"; l.e.. 
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G(x,x') = G 


DIRECT ^ ^TAIL 

+ U 


(4) 


D1R£CT 

where G is nonvanishing on the future light cone of the event x* 

alone, and is given by 


^DIRECT _ ^^^^-1 6[n(x,x’)] . 


(5) 


1/2 

Here A is the "scalarized Van Vleck determinant” (see DeWitt and 

Brehme) , (2(x,x') is the "world-function" of Synge (1960), and 6[(2] is 

the retarded Dirac delta function and is nonzero only when x* lies on 

the past light cone of the event x . (Throughout this paper 6 will 

represent the retarded delta function. The appropriate temporal order of 

the events in will always be x* -< x" -c x where -< n«ans "in the causal 

TAIL 

past of"; cf. Figure 1.) G is a nonlocal term which arises from back- 
scattering of the direct field by the curvature of spacetime. The neces- 
sity to construct approximate Green’s functions is chiefly a consequence 
of the complexity of the tall term. 

For physical situations in which the gravitational field is weak, it 
is possible to choose coordinate systems in which the metric can be wrlt- 
^yv \v “ diag(-l, 1,1,1) and |h^^| « 1 . Under 

such circumstances, to first order, the Green's function can be written 


where 


j^G(x,x') = qG(x,x') + AG , 


qG(x,x') - (4tt) ^ 6[Qfi(x,x')] i (47T) ^6 , 


(6a) 


(6b) 


qQ(x,x’) » ~ )(X^- ^ ) 

= y n 

2 yv 0 0 


(6c) 
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and where AG is the lowest-order contribution to the flat-spacetiiM 
Green's function qG(x,x’) due to the nonzero values of 

Throughout this paper the symbol 6 , standing by itself, will 
denote the retarded Dirac delta function of j|^(x,x'), and similarly for 
6 and 6 : 


6 = 6(qJ2) , 

Z = 6(q?T) , 6 = 6(qQ) , 

(7a) 


, qIT = qQ(x,x") , = qJ^Cx ,x*) • 

(7b) 


Note that 6,6, and 6 are essentially flat-space propagators, while 
6^, and 6^ are 4-dimensional Dirac delta functions of position (eq. 
[3])* A prime on a 6 , 6 , or 6 will always mean derivative with res- 
pect to its argument 

6' E 36/3qJ^ , '6' = 3^/3q?T , 6’ = 36/3qQ . (7c) 

The spatial gradients of are flat-space vectors connecting 

x’,x", and x ; we shall denote them by capital X’s: 








) 


I 







II 


(x^- 


It 

) 





(x^ - x^ ) 


Note that 


5 

= 6'X , 6 , 

* 5’X,, , 

.y 

y ’ ,y 

y 

0^ 

1 „yVy „ 

“ 2 ^ • 

0^ " j n 




y" V' 


(7d) 


(7e) 

(7f) 
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Figure 1 may help one to remen4>er the above conventions. 


b) j^G(x.x') of DeWltt and DeWltt 

DeWltt and DeWltt (1964) derive a weak-field Green's function for 
the vector wave equation of curved—space electrodynamics. Here we sketch 
the obvious specialization of their derivation to the scalar wave equa- 
tion. 

The DeWitt-DeWltt analysis makes use of gauge invariance techniques 
developed by Schwinger (1951) to deal with Green's functions in quantum 
electrodynamics. Specifically, DD take G(x,x’) to be the matrix element 
of 8X1 abstract operator G in a fictitious Hilbert space: 


G(x,x') ■ <x|Gjx'> , (8) 

and they consider the defining differential equation (2) for G to be 
a matrix element of the operator equation 

FG « -1 , with F = -p^(-g)^^^g^’^ f (9) 

where the 1 is the identity operator in the Hilbert space, and the p^ 
are Hermitian operators characterized by the commutation relations 


. [p„.p^]- 0 • <W) 

Taking the variation of equation (9) with respect to the n»tric tensor, 
one finds that: 


AG - G AF G 


where 


, 1 , vl/2r uo VT. yi va yv at-i . 

AF » y [g g + g g - g g J ^gg^Pv 


( 11 ) 

( 12 ) 
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By substituting equation (12) for AF into equation (11) and evaluating 
the matrix element, one obtains 


AG(x,x') 




X Ag M_u 0 „(x",x')(-g")^^^ dV' . 

U I f V 


(13) 


The Green's function, to first order in h^^, can now be obtained by making 
the substitutions 


gpv , ^.uv ^ and G(x,x") « qG(x,x") 

in equation (13), and inserting the resulting AG into equation (6a). The 
result is: 


iGjjd(x.x') 


6 + (4tt) 


-2 


r^"v" ~ 


d^x" 


(14) 


where Iv is the trace-reversed metric perturbation (a single-point func- 
tion; not a bi-tensor) 




= h 


u"v" 




h" 


(15) 


and where 6 and 6 are defined by equations (7). This form of the first- 
order scalar Green's function is stated by DeWitt and DeWitt (1964) without 
proof. 

c) j ^G(x,x') of Thome and KovAcs 

The starting point of the Thorne-KovAcs (1975) development of an ap- 
proximate Green's function Is equation (4). Using the explicit expression 
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for XK Insert equation (4) into the defining equation (2) to 

obtain a differential equation for the tall term: 


yVgTAIL ^ -.(4 tt)"^ 6(0) . 

;jiV ' ' 6 ;>iV 


(16) 


Inverting equation (16) and combining with they obtain an alter- 

native expression for the exact Green's function: 


G(x*x') - (4 tt)”^A^^^ 6(0) + (47r)"^ 


~l/2 ;y" 5^)G(x.x")(-g")^^Vx". 

9 M 

(17) 


It is now possible to make a power series expansion of this equation to 
first order in h^^ . To this end, an expansion for the world function 
0(x,x') can be obtained by approximating the geodesic between the points 
x and X* by the "straight line" 


C-oM-- 


s'*- 


+ XX' 


with 0 < X < 1 


(18) 


The errors Introduced by such an approximation are of second order and 
thus do not affect first-order results. One finds ( without imposing the 
de Donder gauge condition, or any other gauge condition, on 

^ TK Imposed the de Donder gauge condition, but the result is the same 
without it. 


n(x,x') 

Y(x,x') 


-Q(x,x') + y(x,x') 


_ 1 a B 

= T X 


h rt dX 
cxB 
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6(x,x*) 



1 + 2a(x,x') , 


(19a) 

(19b) 


(19c) 


7 


a(x.x’) 


(19d) 


- f I X(l-X) dX . 

<3b 

In equation (19d) Is the Ricci tensor, accurate to first order in 
. Inserting equations (19) into equation (17) one obtains 

lG^jj(x,x») - (4Tr)"^(l + a) 6(^12) + (47r)"^ y<S ' (^Q) 

+ (4tt)“^ f a ,y' 6(-?2) 6ijl) dV . (20) 

; »y u 0 

This expression for j^G(x,x*) actually differs from that of Thorne 
and Rcv&cs in two ways: First, to simplify computations when applying 

the formalism, TK have used a different but equivalent version of the 
tall term (the Integral in equation [20]). However, in their Appendix C 
they prove chat their different version is equivalent Co Che one given 
above. 

Second, TK perform a renormalization (truncation) on the bl-scalar 
Y so as to make their Green’s function and gravitational radiation 
formulas valid for field points very far away from the source. (The method 
of DD cannot reveal the breakdown in the formulas at large distances.) To 
effect a comparison of the two Green's functions we shall deal with 
points X and x* near enough to each other so that both Green's func- 
tions are valid, and that, consequently, the TK truncation is not needed; 
cf. Appendix C of TK. 

It should be pointed out Chat the Green's functions of DD and TK have 
both been developed in a gauge-invariant manner (cf. footnote 1), Of the 
two representations, that of DD has the advantage of formal simplicity. 
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that of TK is superior heuristlcally in that it is a sum of four 
teru , each of which is easily understood physically. Also, it appears to 
us (cf. Kov&cs and Thorne 1976) that the DD form is superior for proving 
theorems, but the TK form is superior for practical cc^utations Involving 
bodies separated by distances large compared to their size. 

III. EQUIVALENCE OF AND 

Assiasing, of course, that there are no errors in the derivations 
sketched above, then in their common domain of validity and must 

be equivalent. In this section we shall present a direct proof of their 
equivalence. This is io^ortant beca\tfe it demonstrates explicitly the rels'^ 
tionship between and and thereby helps explain why the two 

formalisms are powerful for very different aspects of breasst rah lung calcu- 
lations (Kov&cs and Thome 1976). 

By inserting into expression ^14) for^G^^ the identity (see Appendix) 

6 , ,. h ^""'"6 - 5 ^ , 6 + 4 w (^ 6 ’ + ^ 6 ) 

»U f\i 

+ (y6’ + ^^6) - + S5) .,] . (21) 

pV *V 

by using the divergence theorem to convert the volv^ integral of the second 

line into a surface integral, and by perfozming the integration over x" in the 

—4 4 

second term which involves 6 =6 (x,x")» one obtains: 

lGjjjj(x,x') - j^G^j^(x,x') + (4n)”^ | {Jh^ ^ 

+ 5 ,,,.[y6' + 5[y6' A." • <22> 

p V f V M 

Ttie three-surface over which the integral is evaluated is at spatial, null. 
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and Cenporal inflnicy* The surface Integral can be broken into five 
separate terms. Each tern contains the delta function S * or one of 
Its derivatives, which are nonzero except on the past light-cone of the 
event x . Similarly, each of the five terns contains the delta function 
6 , or one of its derivatives, which are nonzero only on the future light 
cone of the event x' . Therefore the integrands are zero everywhere 
except on the intersection of the two light cones. In particular, the 
integr^ds are zero at infinity. Thus the surface integral term in equa- 
tion (22) vanishes md we obtain 

j^Cpp(x,x') • j^G^j^(x,x*) . (23) 


IV. ALTERNATIVE VERSIONS OF THE POST-LINEAR GRAVITATIONAL- 

WAVE-GENERATION FORMALISM 


We now present several alternative representations of -he TK formal- 
ism for calculating the gravitational radiation emitted by fast-mo':lon, 
self -gravitating, veak-field sources. 

In the TK formalism the gravitational radiation is the tlme-dependenc, 
transverse-traceless part of a metric perturbation % which is calculated 
at field points x in the radiation zone with accuracy of second order in 
the internal gravitational field of the source. This second-ord'.r field 
is expressed as a retarded integral over the source region (points x'): 


,h^^(x) - 16 tt 


rU'v’^ .W'V’ 


[(l-h')T^ + (16tt)"^ h'^'^\ptliG(x.x')d^x’ 


( 24 ) 


U V U V — uv 

Here T , t , and h are the nongravitational stress-energy 
LL 

tensor (accurate to second order), the Landau-Llfshitz pseudotensor 
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(accurate to aecond order), and a trace-reversed netric perturbation (ac- 

2 

curate to first order) which satisfy the following coupled equations 

^In TK is denoted , t^ is denoted ^t^ , and h^^ is denoted 

We drop the prefixes to simplify the notation. 


..yv 


^av _ pV ^yo 


uv 


ov 


08 2^^ a, 8 ^ 8,a ‘^o8 ^ ’ 


t. 


of 


(167T) 


-1 {i n“«n. „ 5^ 


X\i 


- (n“^n h®'’ h*^ > + n®\ h“" „ 

yv ,p ,X 'yv ,p tX' 




(25a) 

(25b) 


(25c) 




,o8 


-16tt T 




(25d) 


In their version of the formalism TK insert into equation (24) thel 
own first-order Green's function (eq. 20), with their modified tail 

term. If instead we Insert as given by equation (14), we obtain an 

alternative representation for the second-order field: 


,h (x) - 4 


(U - h') T*''"' . * (16.) -1 


i f[ 5 h''"S" 6 d‘x- dV . 

^ j J ,o ,p 


(26) 
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Equation (26) is, within the constraints of the weak-field approxi- 
mation, completely equivalent to the 2 ^^^^ developed by TK (their eq. [58]) 
so long as the field point x is near enough to the source that gravlta- 
tlonal-redshlft-lnduced phase shifts can be Ignored (see the discussion of 
the truncation of y in §11. c of this paper, and see §1V. c.li of TK) . 
Throu^out this paper we restrict ourselves to such field points. 

The first term of expression (26) is exactly the same as the *'direct” 
plus "whump" fields of the TK formalism. It is the second term which dis- 
tinguishes this formula from theirs. The second term is equivalent to the 
sum of their "focusing”, "transition” and "tail" terms. Since this repre- 
sentation appears to be substantially simpler than that of TK — at least in 
mathematical form — it should be useful in the proof of theorems and in the 
computational details of seme applications. 

Let us now consider the second term of equation (26), which we desig- 
nate 1^^ , to see if it can be manipulated into a more tractable form. 
Reordering the integrations we obtain: 



If ..ow we use the inverted form of equation (25d), i.e.. 


=4 


■ T*'’'’’ 6 d''x' . 


we find 




6 d‘x- 

*Ct ,p 


(28) 


(29) 
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While this Is a quite workable form. It is possible to obtain two 
additional representations by integrating by parts. Considering the 6' 
differentiation first and making use of the de Donder gauge condition, 

“U ^ 

h I ■ 0 (which follows to first order from eqs. [25a, d]), we find: 







h“'^' jjy’v* 



d^x* 


(30) 


The first term, which we identify as A^'^(x), by use of the divergence 
theorem can be converted to a surface integral at spatial, null, and tem- 
poral infinity: 


A^^(x) 




(31) 


The integrand is nonzero only at the intersection of past null infinity, 

. 3 

c9 with the past light cone of x . At that intersection 5 ^,d 

|x'|”^ • |x'|^ • 3/3x°^ , while h^ ^ ]x*| ^ . Thus, whatever may 

be the time dependence of h ^ at c9 , the surface integral vanishes at 
least as fast as |x'|”^ . Consequently, equation (30) becomes 





16h“’e’ 
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where we have made v • of the identity 5 • A nicer expression 

,a p ,cxp 

will result if we interchange the order of differentiation and integration. 
However, in doing so we produce a divergent integral since at large radii 
Ix'l. <Sd^x’ 'v |x'l^ d|x’| , and h®^ ® 'v |x' To avoid the divergence 
we must confine the integral to a finite 4-volume If ^ surrounding the source, 
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then Interchange integration and differentiation, then take the limit as 
the boundary of goes to "infinity" (i.e. , as ^ covers all of space- 

time) : 


rVV 


lira 
3^^ ->00 


4itlJ 




V,. 


J.0(6 


(32) 


In a similar manner, integrating by parts with respect to a' in equa- 
tion (29) and invoking the deDonder gauge condition we obtain 




- ^ I 6 h- - h- - d-x 


,a’6’ 


(33) 


where 




'{6h h 


^x’ 


1 _ 

4it 


. ^’ 3 ' rti’v' , 3 > 

6 h h d . 


(34) 


Again, the surface integral can be seen to vanish at least as fast as 
Ix'l"^ . Thus 


l'^''(x) . 


4tt 


» rcx’3’ ry’v’ ,4 , 

6 h h ^^,3, d X 


(35) 


In summary, the second-order gravitational field, 2^^'^(^)» can be 
written 


,h^^(x) = 4 


I [d-h') h''’“\p,]6d\' 


+ l“''(x) 


( 36 ) 


where I^'^(x) can be represented in the following ways 
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(37a) 


- 

11 

y h 6 


n,u'V 

T 

d^x' dV' 

■1 

(47T)“^ 1 6 

ra'B* 
,a' ^ 


•v' 

d^x' 

- 

{- (47T)“^ 

I 6 H“'e' 

h“' 

'V' , 

»3 


- 

- (4it)"^ 

4 

6 

\a'6’ 


' d«x' 

S3 


'3' 

jjy'V 

'*'**'). J 


dir, -^L\ 
A 

% 





tViv, V // f X rct’3* ry*v* ,4 , 

(x) = - (47 t) j 6 h d x' 


(37b) 

(37c) 

(37d) 

(37e) 

(38) 


All Integrations except that In equation (37e) extend over all of spacetime; 
in (37e) the differentiation must be performed before extending *\r ^ to all 
of spacetime. 

We have given expression (38) for an equation nvimber of its own 

for two reasons: (i) it is badly behaved in point-mass bremsst rah lung calcu- 
lations (see §V below), and (il) it is most naturally thought of, not as 
arising from the DD Green's function (which is the way we derived it), but 
rather from the flat-space Green's function. 

The second point (flat-space Green's function as origin of ex- 
pression [38]) can be seen as follows; Expand the Einstein field equations 
to second-order in the metric perturbation in the manner of TK, and impose 
the de Bonder gauge condition to obtain 






(39) 


(cf. eqs. [9f] and [10b] of TK) . This expression leduces to the first-order 
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field equations (25d) if one neglects the quadratic terras on the right 
hand side. The quadratic terms provide the next- higher-order , nonlinear 
correction to the gravitational field. Equation (39) is written in terms 
of a flat -space wave operator, whereas in TK the same equation was written 
in terms of the wave operator for weakly curved space (TK eq. [10b]). If 
we invert it using the flat-space Green's function ^(x,x') ■ (4tt) ^6 
(rather than using the weakly curved Green's function or iSd^ we 

obtain expressions (36), (38) for • 

We now see that there are at least three distinct ways by which one 
can obtain expressions for the second-order field, corresponding to three 
distinct Green's functions: j^G^jri ^d qG . When one deals with ex- 

tended sources with truly weak fields in their Interiors, the three methods 
and their resulting formulas are all equivalent (aside from the delicate 
issue of validity for field points far far from the source). The DD formulas 
developed here (eqs. [36] and [37]) and the TK formulas (their eqs. [58] and 
[59]), however, have an advantage over the flat-space formulas (eqs. [36] 
and [38]) in their ability to approximate widely separated bodies as point 
masses . 


V. THE POINT-MASS APPROXIMATION TO THE BREMSSTRAHLUNG PROBLEM 

The most important application of the post-linear formalism is the cal- 
culation of the gravitational radiation from a collection of gravitation- 
ally interacting stars — the gravitational bremsstrahlung problem. To this 
problem we now turn our attention. 

It is reasonable to expect that in stellar encounters with Impact 
parameters large compared to stellar radii, the monopole fields of the 
stars should produce the major contribution to the time-dependent part of 
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the second-order field 2 ^^^ from the source, and thence to the gravi- 
tational radiation. Moreover, In calculating the effects of the monopole 
fields, It Is tempting to idealize the stars as point masses. Of course, 
such an idealization violates the weak-fleld assumption of post-linear theory 
since the field of a point mass diverges at the position of the mass. But 
that violation is not important. The important question (formulated so as 
to mesh with the following analysis) is this: 

Consider a near encounter between two stars A and B in which (1) to 
avoid issues of gravitational waves from stellar pulsations, the stars are 
assumed to not pulsate at all; (ii) the stars have weak internal gravity: 

■"a “a • 'b “b 

where r. is the radius of star J as measured in its own rest frame and 

vl 

m^ is its mass; and (iii) the radii of the stars are small compared to 

J 

their Lorentz-contracted distance of closest approach: 

rA+ rg « b/Y , Y = (1- . (40b) 

Here v is the relative velocity of the stars at their point of closest 
approach. Question: Will a monopole, point-mass calculation with the 

formulas of po st-linear theory give (very nearly) the same result for the 
gravitational radia tion emitted as one would get from the same post-linear 
formulas, treating the stars correctly as finite bodies with realistic 
multipole structures? 

Kov&cs and Thome (private ccnnmunication) have proved that the answer 
is "y<=‘s" for their post-linear formulas. In fact, the demand for a "yes" 
answer was a guiding principle in the original development of their formal- 
ism. In this section of the paper we shall show that the answer to the 
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above question Is also "yes" for the DD forms of the post-linear equations 
(eqs. [36] and [37]) t but "no" for the flat-space forms (eqs. [36] and 
[38]). 

Begin with the DD equations. We presume that the equations of motion 
for the source (eqs. [25a,b,d]) are solved In such a manner that the monopole» 
point-mass calculation gives essentially the same stellar trajectories through 
the flat background spacetime as the finite-body calculation. We must then 
check whether expressions (36) and (37) for sensitive to the mono- 

pole, point -mass Idealization. 

The fields in the regions exterior to each of the stars can be repre- 
sented by multipole expansions. Near the surface of star A its monopole com- 
ponent Is greater than or of the order of all other components; at a dis- 
tance r’ A*s monopole component is larger than all others by a factor 
2 

(r'/r^) . Thus, there exists a certain radius R^, measured from the center 
of star A in the rest frame of A, beyond which the field of A can accu- 
rately be considered a pure monopole field. Moreover, we can choose R. such 

A 

that 

r^ « R^ « b/Y . (^Oc) 

A similar situation holds for star B . In the volume integrals (36), (37) 
for the only regions that can be sensitive to a monopole, point-mass 

approximation are the neighborhoods 7|(A) and 7?(B) of stars A and B with 
neighborhood radii and Rg. 

Because the stars always remain physically separated, and because the 
equations governing the first-order field are linear, it is possible 

to split T^'^ and h*^^ into independent contributions from each of the stars 
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(41) 


tMV - - hf + hf . 

VJhen the split (41) is inserted into expressions (36) and (37) for » 

three types of terms result: (i) the linearized field 

jE;;'' - 4| +T»’'’’)6d‘x- , (42) 

(ii) "self-energy terms" which involve h^h^ or hghg or 

, and (iii) "interaction terms" which Involve h^hg or h^Tg or 

« 

The linearized field (42) is obviously insensitive to a monopole, 
point-mass idealization since the field point x lies in the radiation 
zone, which is far outside 7? (A) and 7J(B) . 

When one ignores stellar pulsations (in keeping with our assumptions), 
the self-energy terms lead to a simple renormalization of the active gravi- 
tational mass of each star; they are not time-dependent at infinity, and 
they do not contribute to the gravitational radiation (transverse-traceless, 
time— dependent part of • Therefore, we can drop them from our calcu- 

lation and restrict attention to the interaction terms and the linearized 
field. 

Of the interaction terms, those Involving h^T^ or hgT^ (eqs . [36] 
and [37a]) are the easiest to analyze. The (37a) hT terms can be converted 
into the "hh" form (37b) by the manipulations (27)-(29), equally well in a 
monopole point-mass calculation or in an extended-body calculation. Since 
we will treat all hh terms below, we need not consider the (37a) hT terms. 
The h T term in (36) is best analyzed in the rest frame of star A, where, 

i5 A 

using 
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(43) 


^ -1 
- r 

we can bring it into the form 


6 = 6( fi) K r""*" 6(t' “ t +r - n • x') , r = |x) , n = x/r 




-1 


h'l d^x' 

^B^t'- t-r + n.x' 


(44) 


The integrand is nonzero only inside star A , which is far enough from 
star B that only the nK>nopole field of B contributes significantly to 
h’ . Also, because star A is small compared to the impact parameter 
(r « b), and because the field of B inside A changes on a timescale 

A 

> b/(VY) » r , we can ignore the spatial variation of h' across star A, 

A ® 

and we can also Ignore its time retardation from point to point across A ; 
i.e., we can write 


riiv , -1 ,-*K)no, I ft ^ \ 

2^11 “ ^B t’ « t-r) 


t'?'''' a - 

A 


= -4r"^ = 0, t’ - t-r) 6^ 6^ . 


(45) 


This is precisely the same result as one would get from the monopole, point- 
mass approximation. The interaction term involving h^T^ in expression (36) , 
when analyzed in the rest frame of star B , gives a similar result. Thus, 
all the hT interaction terms are amenable to the monopolc, point-mass ap- 
proximation. 

There are many interaction terms involving A large number come 

from the t^ ^ term of expression (36) (cf. eq.[25c]); two come from the term 

II ’v) ' 

following t^ in (36); and two each come from expressions (37b,c,d,e). 

iaij 

Each such term involves an integral over all spacetime (all x*). The only 
portions of the integrals which could possibly be sensitive to the monopole, 
point-mass idealization are the portions which come from the non-monopole 
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regions 7KA) and 7?(B); and because our analysis Is Insensitive to the 
change of names A B , It Is adequate for us to consider contributions 
from 7? (A). 

Each of the extended*~body h^hg Integrals over 7^ (A) can be split Into 
two parts: the contribution from the Interior of star A, 7{^(A), and the 
contribution from its exterior, 7|j,(A): 


7?j(A): 0 < r' < r^ , 7?g(A); r^ < r’ < 


u’v* 

Because we ignore pulsations of A , Is Independent of 

— ct ' 8 * “0 * 0 * 

first order the only nonzero component of h^ is h^ 

-0 ' 0 ' 

terior of A , h. and its spatial derivatives have magnitude 

A 


(46) 

t ' , and to 
In the in- 


2 — fi'n* 

\ '’a ‘'a ,J'k' ’JitA). (47 a) 

-O'O' 

In the exterior of A we can expand in multipoles, obtaining 


00 

hj ^ * (4m /r') {l + ^ a»(0',4)')(r /r’)*^} in 7?„(A) , 


(47b) 


i=2 


where the are all of order unity. Because b >> r^ , the field of B 

inside 7?(A) is (very nearly) equal to B's monopole field at the center of 
A, plus a fractional correction of order yr'/b: 



I 


h^'“' 

B mono 


(t',x'= 0)(1 + 0(yr Vb)} . 


(48a) 


The motion of star B past star A causes the field of B near A to vary 

— u * V * 

on a timescale At' > b/vy >> R. , leading to a time derivative of h_ 

A a 

given by 
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(48b) 


where t'* 0 Is the momait of closest approach » near which time 


0/3t')hj;’^'(t’,x’- 0) 

mono 


I7TV l^’l 

0/3t’) b/n.x’- 0) ** 


vy 


(48c) 


The error terms In (48b) are much smaller than the leading term at all times 
except jt’l r'/vy; and the error terms for |t*| < r'/v are negligible 
compared to the leading term for nearby times |t'| b/vy . Because the 
field point x is very far outside the source region 71(A), we can write 
the propagator 5 in the form (43); uid we can write 


.j' 




,j'k' “ 


(49) 


Equations (43) and (46)-(49) are the foundation for evaluating the h.h„ 

A B 

interaction integrals. By inserting them into each interaction integral in 
turn, one can verify that all the interaction integrals in (36) and (37) are 
insensitive to the point-mass idealization. Consider, for example, the 
Interaction integral 


J^'^“(X) E 


. ra'6’ rM*v' 

^ ^B ,6' 


d x' 


(50) 


71(A) 

which comes from expression (37c). By inserting expression (43) for 6 and 

—Cl* 6 * 

integrating over t' , and by using the fact that h^ is independent of 
t' and is zero unless a' *3' « 0, we bring this into the form 


J"^^“(x) = 


-l.a 


71(A) 


rO’O' rU’V , 

"a ^”b ,0'^t'-t-r+n 


d\' 


(51) 
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We next split the integral into interior and exterior contributions; ve 
use (47a) and (48b, c) in the interior, «id (47b) and (48b, ^ in t'le exterior, 
thereby obtaining 


tWVO 


(x) 


ISTTm 


A 2 3 cW’v’ 

'a “ 


•K— h^ '' (t*« 

^ "b mono^*" 


t-r,x 


0 ) 


1 + 




, (52a) 



ISTTm. 


2 3 rW’v’ 


Vir's 


mono 


(t'= t-r, X*- 0) 


X 



(52b) 


Notice that (i) the exterior contribution dominates over the interior 

contribution by a factor (ii) except for the "error .erms" 

the exterior contribution is precisely the result which one would obtain 
from a monopole, point-mass calculation; (ill) the error terms are negli- 
gible at all times except |t-r| < R^/v when the dominant, "point-mass con- 
tribution" is going through zero; ( iv) even when they dominate (for 
|t-r| < R^/v) the error terms are negligible, 0(yR^/b), compared to the 
point -mass contribution at nearby times (|t-r| b/vy). These facts show 

that, for all practical purposes, the Interaction integral (50) is Insen- 

3 

sitive to a monopole, point-mass idealization. 


1 3 

One can show that because of their angular dependences Inside the d x* 

integral, the error terms shown explicitly in (52) actually vanish. Thus 

the dominant errors are of even smaller order than indicated in the 0[ ] 

expressions of (52). 


All other h.h_ terms in expressions (36) and (37) can be handled 
A D 

similarly, giving the same conclusion. 
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Thu8| for the DD version of the post-llneer forBallsn (eqs. [36] 
and [37]) the answer to our question (italicized sentence following eq. 
[40b]) is "yaa". The aonopole, point-iBaso approxlnatlon is valid. 

Not so for the flat-space version of the post Inear fomalisn 
(eqs* [36] md [38]): Tlte second derivatives which occur in expression 

(38) wreak havoc with the nxmopole, point-aiass Idealization* To see this, 
consider the following interaction integral, which coues from expression 
(38): 




7](A) 


6h 


uV’v’ 


t'3’ 


A .o' S’ "b 


d^x' 


(53) 


Evaluate the Integral over t' in the rest frame of A , using expression 

“U'V* 

(43) for 6 and using the fact that h^ is independent of t and is 
zero at first order unless y' - v' - 0’ ; the result is 


. r' 


■1 xU xV f cO'O’ rJlj'k'* .3 , 

^0 S J ^A ^t’“t-r + n*x' ^ 


Tlie contribution to from the interior of star A is obtained by in- 

serting (47a) and '48a) into (54), and Integrating over r’ r^ ; the 
result is: 

(16tto /r)6}J hj ^ (t'- t - r, x'« 0) • [1 + 0(Yr*/b) ] . (55) 

X A u u B mono - a 


This in. frnal contribution does not go to zero as the size of star A , r^ , 
is shrunk to zero. Independently of r^ , it is conq>arable in mapiitude 
to the total gravitational-wave amplitude — and for observers not in the 
rest frame of A it Is an indispensible, non-negligible contributor to 
the gravitational waves. A monopole, point— mass calculation will miss this 
contribution. 
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VI. CONaUSION 


Th«re are three different Green's functions which one can use to 
calculste the second-ordcr gravitational field of post-linear theory: 
1®TK* l^DD* 0^’ first-order Green's functioi's j^G^ and are 

coiq)letely equivalent. ^G^ leads to the TK fomulas for 2 ^^^ (their eqs. 
158] and (591), while j^G^p leads to the formulas of this paper (eqs. [36] 
and [37]). Both sets of formulas (TK «id DD) remain valid if one idealizes 
the stars of a bremsst rah lung calculation as monopole point masses. However, 
the formulas for 2 ^^ which are obtained from the zero-order Green's func- 
tion qG (eqs. [36] and [38]), although valid for extended bodies, are not 
valid in the monopole, point-mass idealization of a bremsst rah lung calcula- 
tion. 
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APPENDIX 


The proof of the equivalence of the Green's f\mctlons of DD and IK 
Is based upon the identity (eq. [21]): 

T ® v" " ^ +4 it(y 6' + a6) 5^ 

+ [6h^"^"6 .,+ ,.(y 6' +a6) -ri^^6(y6'+ a6) ,,] . (Al) 

|V fM 

(See eqs. [1] and [7] and fig. 1 for notation.) We shall now prove this 
identity. 

If we carry out the differentiations of the rotmd and sqtiare brackets, 
make use of the defining differential equation for the flat-spacetime Green's 
function 

6 „ = 2 6" + 46' = -4ttS^ (A2) 

,y V 0 

and the fact that 

a 6 = -4ir a 6^ = 0 

which follows from equation (19d) since a(x',x') = 0, the proof of the 
identity is reduced to showing that 


K = X^„ h^ ^ 6' + h^ ^ \y"v" " \y"^ 


■2V>' Xj^„6" - y n“'' - 2S’'' X|J,. 6' - 0 . 


Using the definitions of a and Y given by equations (19b, d), it is 
relatively straightforward to show that 
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u" _ o:.y" ' 




,11,. II 

y" .V 


(A4) 


and 


*■ It*' — ~ 1 "“u" ~V)" 

Y-'' X^„ -Y + ^X'^ X'' h^„^„ 


(AS) 


Making use of pquations (A4) and (A5) in (A3) we find 


K - - h"6’ - X^„X^„ h^"^" 6" - Y(26"+n^^ 


If we now employ the identity 


6 „ = X „X „ 6" + n 6' 

,y"v y V yv 


(A7) 


as well as the definition h^ ^ = h'^ j h" and equation (A2) , we 

find 


K = 2tt h" 6^ - ^(26" + ^\y"v"> • 


(A8) 


Using equation (A5) for y and the identity 


6’ „„,,II = 2 Jl6"' + 46" 


,y V 


0 


(A9) 


we can rewrite (A8) as 


K 


= 2tt h" 6^- (y’^"- y X^.i h^”^") X ,„(65" + 2 „f26"’ ) . 


(AlO) 


By differentiating (A2)we obtain the identity 


X^„(66" + 2 Qfi6"') = -4 tt6^ . 


(All) 


Making use of this in (AlO), we obtain 
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K 


2tt h" 6^ + j X^„ h 

27t h" 6^ + “ X^m 

- 47r(?^"y.. - I h" - i X^„ h 


y"v” 

hU-v 

y"v" 




(A12) 


By virtue of the properties of the Dirac delta function and the fact that 

H 

*^(x*,x') = 0 , X^„(x',x') = 0 , and y’’^ ^.i(x',x') = h" , the second 
term vanishes and the third term reduces to -27Th” 6^ which cancels the 
first term. Thus K = 0 and the identity (Al) is established. 


28 



REFERENCES 


DeWitt, B. S 
DeWltt» C. M 
Kov&cs, S. J 
Schwinger, J 
Sjmge, J. L. 
Thome, K. S 
as "TK" 


. and Brehme, R. W. , Ann . Phys . (USA) , 9, 220, 

. and DeWitt, B. S. 1964, Physics , 1, 1; cited In text as "DD**, 

. Jr. and Thome, K. S, 1976, Ap , J . , paper in preparation. 

. 1951, Phys . Rev . , 82, 664. 

1960, Relativity: The General Theory (Amsterdam: North Holland). 
. and Kovacs, S. J., Jr. 1975, Ap . J . . 200, 245; cited in text 


29 



FIGURE CAPTIONS 


Figure 1. A mnemonic diagram for helping one to remember: (i) the tem- 
poral order of the events x’,x",x; (ii) our notational conventions 
for 2-point functions (eq. [1]) and especially for qQ, qQ, qA, 5 » <S, 
6 (eqs. [7a, b]); (iii) our definitions of X , X , and X (eq. 
[7d]); and our formulas for the spatial gradients of the propagators 
6 , 6 , 6 (eq. [7e]). 
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